g-Factor anisotropy of hole quantum wires induced by the Rashba interaction 
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We present calculations of the g factors for the lower conductance steps of 3D hole quantum 
wires. Our results prove that the anisotropy with magnetic field orientation, relative to the wire, 
originates in the Rashba spin-orbit coupling. We also analyze the relevance of the deformation, 
as the wire evolves from 3D towards a flat 2D geometry. For high enough wire deformations, the 
perpendicular g factors are greatly quenched by the Rashba interaction. On the contrary, parallel g 
factors are rather insensistive to the Rashba interaction, resulting in a high g factor anisotropy. For 
low deformations we find a more irregular behavior which hints at a sample dependent scenario. 

PACS numbers: 71.70.Ej, 72.25.Dc, 73.63.Nm 



I. INTRODUCTION 

Spin-orbit interactions in semiconductor materi- 
als offer interesting possibilities of spin control in 
nanostructures.- Among them, the Rashba interaction 
that originates in externally applied electric fields is most 
promising due to its tunability. In this work we prove 
that the Rashba interaction is an important source of 
spin anisotropy in hole quantum wires. This anisotropy 
manifests in large differences between the energy split- 
tings for magnetic fields parallel and perpendicular to 
the wire. Our calculations show that in the presence of 
Rashba interaction the perpendicular field becomes much 
less effective in generating spin splittings than the paral- 
lel one. This effect is favored by the deformation of the 
quantum wire, i.e., anisotropy increases when the wire 
evolves from 3D towards a more flat quasi 2D geometry. 

In semiconductor hole systems like p-type GaAs nanos- 
tructures transport is mediated by holes in the valence 
bands. As compared to electrons, holes are characterized 
by a spin 3/2, besides a sign difference in charge. The 
corresponding fourfold discrete space is a source of qual- 
itative differences with respect to the more usual twofold 
spin of electrons. In 2D hole gases different splittings for 
normal and in-plane fields have been observed, as well as 
for different in-plane orientations.- By further confining 
the hole gas it is possible to generate nanostructures with 
the shape of quantum wires. In this case, the splitting 
varies, in principle, with both wire and magnetic field 



orientations 
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There are few theoretical analysis of the spin splittings 
in hole quantum wires.^^- Although the Rashba interac- 
tion was usually not taken into account, this situation 
changed in some recent worksiiSiii Indeed, Quay et alw^ 
have observed the formation of a spin-orbit gap induced 
by the combined action of magnetic field and Rashba 
coupling in a hole quantum wire, while Chesi et aL— 
have studied, both experimentally and theoretically, the 
spin resolved transmission of a quantum point contact 
fabricated in a 2D hole gas. In the latter, the Rashba 
interaction is shown to favor a band crossing at finite 



wavenumber that can be manipulated with an external 
magnetic field. In agreement with our results, this cross- 
ing is obtained in a multiband description of the hole 
states. It can also be explained within a restricted single 
band description adding a cubic Rashba term. 

In this work we have focussed our attention on the 
structure-inversion-asymmetry (Rashba) splitting since 
this is known to be the dominant source of spin-orbit 
coupling in GaAs. The bulk-inversion-asymmetry (Drcs- 
selhaus) is much smaller and has a minimal effect on 
the energy bands. We will show that in a hole quan- 
tum wire oriented along x' the Rashba interaction due 
to asymmetry in the growth direction (z') causes a large 
difference between parallel {x') and perpendicular [y') g 
factors of the wire, as deduced from the _B-induced split- 
tings of the conductance steps. This anisotropy is due 
to the quenching of the splitting when B is along y' and 
the wire flatness is large. For smaller deformations the 
situation is less clear due to a non monotonous evolution 
of the splittings that may result in a sample-dependent 
scenario. 



II. MODEL 

We describe the anisotropic kinetic energies 
the holes in a 4-band kp model. Introducing a spin dis- 
crete index 77 = 3/2, . . . , —3/2 and following the notation 
of Ref. [l[ the diagonal terms read 
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(71 + c„72)fc^ + (71 - 2c^72)fc' , (1) 



where c±^/2 — 1 and c±i/2 — —1- In Eq. ([T]) 71 and 72 

are the kp parameters, k is the 3D wavenumber and we 
have also defined fc^ ~ k'^ + ky. The nondiagonal kinetic 
terms are 
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where k± — kx ± iky and K — k^ — ky. We only refer to 

contributions in the upper triangle of matrix 7^^^™^ since 
the remaining ones can be inferred from the Hermitian 
character of the matrix. In all calculations discussed be- 
low we have used numerical values for the kp parameters 
7's corresponding to GaAsji 

The wire confinement is represented by a deformed 2D 
harmonic oscillator. Assuming the wire is oriented along 
x' while transverse and growth directions are given by y' 
and z' , respectively, it is 
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The adimensional parameter a of Eq. (|3]), corresponding 
to the ratio of confinement strengths in z' and y', controls 
the flatness or 2D character of the wire. The direct cou- 
pling with the magnetic field B is given by the Zeeman 
term 



-2k^bB ■ J , 



(4) 



where k is a kp parameter, /is represents the Bohr mag- 
neton and J is the angular momentum operator for a spin 
3/2. Finally, the Rashba interaction is described by 



•H^^) = (fc X 7?) • J 



(5) 



where we defined a vector constant TZ = a£, related to 
the effective electric field £ and kp parameter a.A We 
shall treat 7?, as a two-parameter vector with dominant 
component along the growth direction, i.e., TZ = TZziUz' + 
TZy'Uyi with TZz' > TZy' ■ 

In the presence of a magnetic field, the orbital ef- 
fects of the field are taken into account by means of the 
substitution k — > — iV — -^A with the vector potential 
A= {—yBz + zBy, —Bxz/2, B^y /2). In this process, Her- 
miticity is enforced in the cross terms by using the sym- 
metrized forms, such as k^ky — >■ {k^ky + kyk^) /2. Sum- 
marizing all contributions the total Hamiltonian reads 
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The wire Hamiltonian eigenvalues can be labelled with q, 
a real number representing the longitudinal momentum 



and an index 7 = 1,2,.. 



as 
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where ei{q) are the discrete energy bands of the nanos- 
tructure. The eigenvalues are ordered as ei{q) > £2(9) > 
. . . since the spectrum is not bounded from below due to 
the negative kinetic terms. 

We have obtained the solutions of the eigenvalue prob- 
lem given by Eq. ([7]) by discretizing in harmonic oscillator 
states for the two transverse oscillators along y' and z' , 



\Iq) = 1""^^) ' 
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where n, m = 0, 1, . . . represent the number of quanta in 
each oscillator, respectively. The resulting matrix eigen- 
value problem reads 



J2 {n'm'rj'\n{q)\nmrj} C^^^l = ei{q) C\: 



(9) 



In practice the number of oscillator states in expan- 
sion Eq. ([U can be truncated once convergence of the 
results is ensured. The results shown below are well con- 
verged and they have been obtained including the lower 
20 oscillator states in each direction. In Appendix |B] a 
precise discussion on the relevance of the basis truncation 
is given. 



III. RESULTS AND DISCUSSION 

As illustrative examples. Fig. [1] displays the energy 
bands of selected cases. As is well known, the Rashba 
interaction causes a characteristic band structure easily 
recognizable by the pairs of subbands crossing at q = 
and with maxima at opposite q values (left panel) . These 
maxima correspond to band energy minima for the case 
of electrons. In the presence of a magnetic field, when this 
points along the wire (x', central panel), an anticrossing 
of the bands appears at q = 0. This anticrossing may 
lead to anomalous conductance steps, similar to those 
recently measured in Ref. In Fig. [1] this behavior can 
be seen for (£', q) w (— ll?iajoi 0). For B in the transverse 
direction (?/', right panel) the band crossings persist, but 
the two central maxima for each pair of bands are shifted 
differently in energy, the band structure becoming asym- 
metric with respect to q inversion. 



The B-induced modifications of the band structure, as 
seen in Fig. [H cause a change in the conductance of the 
wire. This modification of the conductance, in the limit 
of weak magnetic field, is conveniently summarized by a 
number called the g factor of each conductance split level. 



At = 0, time reversal invariance of the system causes 
the conductance G to increase in steps of 2Go as the 
Fermi energy of the leads is reduced, where Go = /h 
is the conductance quantum. The evolution of the wire 
conductance with energy can be understood if we imag- 
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FIG. 1: Energy bands for a — 64, TZ^' = 2.6hLJo£o and TZyi = 0. Left panel is for B = while center and right ones are 
for /ifl-B = O.lftcJo in the parallel and transverse directions, respectively. The wire is oriented along (—2,3,3) and the growth 
direction is (3, 1, 1). 



ine a horizontal line, indicating the position of the Fermi 
energy, in the left panel of Fig. 1; as this line is moved 
to lower energies it sweeps the band maxima always in 
pairs, each maxima corresponding to an increase of Gq 
in the conductance for hole transport. The result is the 
typical staircase conductance, with step heights of 2Go. 
A similar procedure for the central and right panels of 
Fig. [1] convince ourselves that intermediate half steps in 
conductance are caused by the magnetic field. They are 
smaller than the full steps and proportional to the inten- 
sity of the magnetic field. 

The scenario we have just sketched is explicitly shown 
in Fig. [21 highlighting the conductance half steps at odd 
multiples of Gq. Notice that the energy span varies for 
each specific conductance half step. In the limit of weak 
magnetic fields we can conveniently summarize the B- 
induced A^th half step in the conductance, appearing be- 
tween steps at 2{N—l)Go and 2NGo, in terms of a single 
number called the g factor. As this number depends on 
the conductance step and the magnetic field orientation. 



we use the notation g^^^"^ and g'"f ' to indicate the g factor 
of the A^th step, for B along x' and y' , respectively. Of 
course, other orientations are in principle possible, but 
we will restrict first to these two as they are the rele- 
vant ones in the measurements of spin hole anisotropy. 
In Appendix |^ we will briefly mention the behavior for 
z'-oriented field. 

Our precise definition of the parallel-field g-factor is 
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where is the energy range for the A^th half step in 

a magnetic field B. In Eq. pO)) . the factor 3 in the de- 
nominator is introduced by conventions^ The definition 
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of g^" , for magnetic field along y' , is obtained simply 
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by A^'^ in Eq. (H 



FIG. 2: (Color online) Conductance traces for the band struc- 
ture in parallel and perpendicular magnetic fields of Fig. [T] 
For comparison, the 5 = conductance is displayed as a thin 
line. 



Figure [3] displays the perpendicular (lower row) and 
parallel (upper row) g factors for the lower conductance 
steps, as a function of the wire deformation a and for 
different values of the Rashba coupling TZz' ■ These are 
the main results of our work. They were obtained for a 
specific wire orientation and direction of crystallographic 
growth (z') taken from the experimental works of Dan- 
neau et al^ and Koduvayur et al^ We have checked, how- 
ever, that a qualitatively similar infiuence of the Rashba 
intensity and confinement deformation are obtained as- 
suming other arbitrary orientations. The g factors show 
a general tendency to decrease as a increases, except for 
smaller deformations (a < 100) for which g may increase 
or even show irregular behavior in some cases. Focussing 
first on g\\ , we notice that this component does not change 
significantly when the Rashba intensity increases, spe- 
cially at large a's, for which the results are almost over- 
lapping in the upper panels of Fig. [31 Very remarkably, 
however, for magnetic field in the perpendicular direction 
small variations in TZz' are enough to strongly modify the 
values of This is more clearly seen in Figure [31 which 
displays the dependence with Rashba coupling intensity 
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FIG. 3: (Color online) Parallel and perpendicular g factors as a function of wire deformation for different values of the Rashba 
strength: TZ^i = (circles), 1.5huJo£o (triangles) and 2.6htjJo£o (squares). Upper and lower rows are for parallel and perpendicular 
fields while columns from left to right correspond to increasing conductance half step A'^ (see text). The results for A*' = 1 are 
not shown due to their similarity with the displayed N = 2 case. The orientation of the wire is the same of Fig. [T] 
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of the g factors. 

There is a general Rashba-induced quenching of 
in Figs. |3] and SI quite conspicuous for iV = 4 and 5. 
This effect is so strong that it can reverse the relative 
importance of f/|| and g±; from g± > g^^ when TZ^i = to 
g± « g\\ for increasing TZz' (> 2.5hujo^o, Fig. S]). With 
the chosen values of TZz' we even find a range of a's for 



which essentially vanishes. It is interesting to point 
out that a similar quenching of conductance plateaus in 
transverse field was discussed in Ref. [l^ for parabolic 
wires with electron conduction, as opposed to the present 
hole conduction. In both cases the Rashba spin-orbit 
coupling is the underlying mechanism. 
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Turning to the comparison with experiments, this is 
somewhat complicate due to the sample dependence. In 
general, however, a large g-factor anisotropy between 
parallel and perpendicular orientations has indeed been 
observed in Refs. O-IBl. This was generally attributed 
to a preferential orientation of the spins along the wire 
for strong confinements. Our results prove with de- 
tailed calculations that the Rashba interaction for holes 
is the specific mechanism allowing the appearance of this 
anisotropy. As this interaction is sample dependent and 
may vary with external field, our results also predict that 
the hole g factors may be tunable to a certain degree, 
what may be relevant for spintronic applications. The 
experimental values of wire deformation a are somewhat 
uncertain in general, which is an additional source of dif- 
ficulty for comparison. In general, however, experimental 
wire deformations are a < 100, which in our calculations 
corresponds to a regime with rather large fluctuations 
(Fig. 131). Only for larger a's the value of is consis- 
tently below ^y^-* at high enough TZz'- We believe that 
detailed comparison in this regime is quite involved due 



to the fluctuations. On the other hand, these sharp vari- 
ations of in the small-a regime and of g± at all a's 
can be seen as a manifestation of magnetoconductance 
tunability via the Rashba coupling. 



IV. A TWO-BAND MODEL 

A more transparent physical interpretation, comple- 
menting the above numerical results, can be obtained in 
a simplified model based on only two bands. Focussing on 
the I-th intermediate half step having conductance I Go, 
with / = 1,3,..., we select the two states / and / -I- 1 at 
a given q, {\Iq)o, 1(1 + ^)q}o}, where the zero subscript is 
indicating absence of a magnetic field. These two states 
are the basis in which the effect of the magnetic field in 
different orientations will be described. 

Let us assume that the B-field Hamiltonian may be 
split as 
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where T-l^^^ is the Zeeman energy defined above in Eq. ^ 
and Ho is the zero field Hamiltonian in Eq. Notice 
that Eq. (fTTj) neglects orbital field effects, a simplifying 
assumption motivated by the qualitative nature of the 
present two-band model. 

The zero-field energy bands, given by 



'Ho{q)\Iq)o = e/o(g)|/g)o , 



(12) 



are assumed known; such as those displayed in the left 
panel of Fig. [1] for a specific confinement and Rashba 
intensity. In presence of a magnetic field the modified 
energy bands are the eigenvalues of the matrix 
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FIG. 5: (Color online) Wavenumber dependence of the matrix 
elements entering Eq. (|15|l and (|16|) for 1 = 7. Circles, trian- 
gles and squares are for increasing values of the Rashba in- 
tensity TZ^' = O.lhujoio, 1.5hujo£o and 2.6?icjo^o, respectively. 
(Other parameters: a = 150, B = 0.) 



In a parallel field "H*-^-* oc Jx' and, for this case, 
we have found that the 7jS vanish. This is reminis- 
cent of the behavior of conduction electron wires, where 
the spin textures also show a vanishing integrated spin 
along the wireji^ In parallel orientation the band ex- 
trema are at g = (see Fig. [TJ middle panel) for which 
£/o(0) = e(7+i)o(0) due to Kramers degeneracy. Under 
these conditions we find from the two eigenvalues of the 
matrix in Eq. ()13|) that 



^11 



o(/0|J,H(/ + l)0)o 



(15) 



That is, the parallel g-factor is determined by the transi- 
tion matrix elements of the parallel spin component be- 
tween the Kramers degenerate states at q = 0. The upper 
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FIG. 4: (Color online) Parallel (left) and perpendicular (right) 
g factors as a function of Rashba coupling intensity for dif- 
ferent values of the deformation a. The results correspond to 
the N = 5 conductance half step. 



panel of Fig. [5] shows this transition matrix element for 
1 — 7. Notice that for g « the transition matrix ele- 
ment is not depending on the Rashba intensity, thus ex- 
plaining why the parallel g-factor is not strongly affected 
by the spin-orbit coupling. 



B. Perpendicular field 

For H*-^-' (X Jy' the band maxima are shifted in op- 
posite directions for positive and negative q's (Fig. [T] 
right panel). This implies that the energy difference de- 
termining the g factor corresponds now to states with 
opposite wavenumbers, say qm and — g„i. For nonzero 
qm the two states Sio{qm) and e(/+i)o(9m) are nonde- 
generate and, for a sufficiently small field, we should 
have 6 « £i{qm)ili in Eq. ([T^. As a matter of fact, 
we find that 5 actually vanishes for the perpendicular 
field. This is the regime of non-degenerate first-order 
perturbation theory with modified energies £io{qm) + 7/ 
and £(/+i)o(9m) + 7/-i-i- With the explicit definition 
of the 7's and noting that Sioiqm) = £ioi—qm) and 
7/(9) — for any q (Fig. [S]) the perpendicular 

g factor reads 
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It seems natural that in y' orientation the g factor 
is simply proportional to the expectation value of Jy'. 
Figure [5] shows the variation of this expectation value 
with the wavenumber and the Rashba intensity. Notice 
that typically —0.5 < qm < 0.5, i.e., the maxima are lo- 
cated in the central part of Fig. [5] lower panel. When 
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the Rashba intensity TZz' increases there is a severe re- 
duction of {Jy')o in absolute value for this central region. 
This is the mechanism by which the Rashba interaction 
quenches the transverse g factor; namely, by means of a 
strong reduction of the transverse y' spin component. 

For strong spin-orbit coupling the expectation values of 
all three components of the spin vector at zero magnetic 
field, ( J)o, vanish; a manifestation of the spin randomiza- 
tion induced by the Rashba field TZ. In y' orientation this 
induces a quenching of the g factor through Eq. but, 
quite remarkably, Kramers degeneracy at zero wavenum- 
ber keeps the parallel g factor almost unaffected by virtue 
of the transition matrix elements in Eq. (|15p . 

The g factors obtained from Eqs. and nicely 
agree with the results from the full diagonalization when 
orbital effects of the magnetic field are also neglected in 
the latter. The comparison with the complete model, 
results of Fig. [31 is less good; the trends are qualitatively 
reproduced but differences may be as large as a factor 
two. Orbital effects of the field are thus quite important 
for a precise analysis. 



V. CONCLUSION 

We have attributed the anisotropy of magnetotrans- 
port g factors in hole quantum wires to the Rashba in- 
teraction. When the wire deformation and Rashba inter- 
action are both large enough (a > 100, TZ^' > 2.5huJoio) 
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is greatly quenched by the Rashba interaction and 
is almost unaffected. For lower wire deformations 
(a < 100) we find a fluctuating, sample dependent be- 
havior of the g factors. 
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Appendix A: Field along z' 

Experimental g factors are usually obtained for mag- 
netic fields in the x'y' plane, either in parallel (x') or 
perpendicular (y') direction with respect to the wire. For 
completeness, in this Appendix we discuss in a qualita- 
tive way the effects of the magnetic field when this points 
along the growth direction z' . The energy bands are sim- 
ilar to those of the x' orientation (middle panel of Fig. 
[1]): they are symmetric respect to 5- inversion, with anti- 
crossing points at q = 0, although the i3-induced split- 
ting is much stronger. This enhancement agrees with 
experiments^ and is surely due to the important orbital 
motions induced by the field in this geometry. We thus 
obtain gz' > g\\ , where gz> and g\\ denote the g factors for 
z' and x' fields, respectively. 



Looking at the Rashba-field dependence, gz' behaves 
similarly to g± (along y'): it decreases with increasing 
TZz' but does not vanish for the maximum value we have 
taken (2.6?ia;o^o)- For strong wire deformation the sat- 
uration value corresponds to gz' ~ 5, while for in-plane 
magnetic field it corresponds to ^n^^ w 1.5 (see Fig. [3]). 
For small values of a the behavior of gz' is less regular, as 
for the other orientations, but it tends to increase with 
a. Within the two-band model of Sec. IV we expect 



o(/0|J,.|(^ + l)0)o 



(Al) 



which is equivalent to Eq. (|T5|). replacing J^,' — >■ Jz' , and 
is now depending on the value of the Rashba intensity. 



Appendix B: basis truncation 

This Appendix discusses the relevance of the trunca- 
tion of the number of oscillator states for the y' and z' 
oscillators. It is usually assumed that the confinement 
allows the truncation to the lowest, or few lowest, states. 
Here we explicitly check this quantitatively for selected 
values of a, the ratio of the two confinement strengths. 
We restrict, for simplicity, to the 5 = case with strong 
Rashba coupling in the growth direction. 

Figure [B] displays the evolution of the band structure 
when a) increasing Nyi and Nz' sequentially from left to 
right panels; and b) increasing the deformation degree 
a from top to bottom panels. The right column shows 
results that are very close to physical convergence. Look- 
ing at the successive band crossings at q = 0, we notice 
that the truncation to {Nyi,Nz') = (1,1) grossly overesti- 
mates the energy separation between pairs of bands in all 
cases. It is remarkable that for increasing flatness degree 
the (1,1) truncation deviates more and more from the 
right column. This is a consequence of the intersubband 
couplings induced by the kp and Rashba Hamiltonians: 
at least a few bands in the shallow oscillator {y') are es- 
sential even for large a's. 

More reasonable results are found for {Ny',Nz') = 
(10,1), although the differences with the (10,10) basis 
are still large quantitatively. In this case, however, in- 
creasing a improves the quality of the description since 
intersubband coupling is allowed at least in y' direction. 
Finally, the (10, 2) results are close to the converged ones 
and only the insets reveal that sizeable differences are 
present at intermediate or low values of a. These dif- 
ferences are small in the behavior of the upper bands 
and become more and more important as the energy is 
reduced. From this analysis we conclude that for our 
present purpose, namely the description of magneto- g- 
factors of several successive conductance steps, it is es- 
sential to include enough oscillator bands in both y' and 
z' oscillators. 
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(JVy,JVJ = (l,l) (10,1) (10,2) (10,10) 
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FIG. 6: Evolution of the energy bands for selected numbers {Nyi , N^i) of oscillator states in the matrix discretization (columns) 
and aspect ratios a (rows). The gray colour results are qualitative, indicating that the corresponding energy regions are full 
of bands. The insets in the rightmost columns show the details of those dense band distributions. Parameters: B = Q, 
TZ^i = 2.6?uooio, TZy' = 0, growth direction (001) and wire orientation (110). 
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